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Abstract: We study the spontaneous breakdown of SO (10) symmetry in the IIB matrix 
model, a conjectured nonperturbative definition of type IIB superstring theory in ten di- 
mensions. Our analysis is based on a Gaussian expansion technique, which was originahy 
proposed by Kabat-Lifschytz and applied successfully to the strong coupling dynamics of 
the Matrix Theory. We propose a prescription for including higher order corrections, which 
yields a rapid convergence in a simple example. This prescription is then applied to the 
IIB matrix model up to the third order. We find that the 'self-consistency equations' allow 
various symmetry breaking solutions. Among them, the solution preserving SO (4) symme- 
try is found to have the smallest free energy. The value of the free energy comes closer to 
the analytic formula obtained by Krauth-Nicolai-Staudacher as we increase the order. The 
extent of the space-time in the 4 directions is larger than the remaining 6 directions, and 
the ratio increases with the order. These results provide the first analytical evidence that 
four-dimensional space-time is generated dynamically in the IIB matrix model. 
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1. Introduction 

One of the recent excitements in string theory is the appearance of concrete proposals 
for its nonperturbative definitions. The IIB matrix model [1], which was conjectured to 
describe type IIB superstrings in 10 dimensions, has a particularly simple form. It is a 
supersymmetric matrix model, which can be obtained from the zero-volume limit of lOd 
SU(A^) super Yang-Mills theory. The space-time is represented by 10 bosonic matrices, 
and hence treated dynamically. This allows us in particular to investigate the possibility 
[2] that our 4d space-time is generated dynamically in superstring theory in lOd. Since the 
model has manifest SO(IO) symmetry, the emergence of 4d space-time requires the SO(IO) 
symmetry to be spontaneously broken. However, the absence of SSB has been concluded 
in various simplified versions of the IIB matrix model [3, 4, 5]. (See also Ref. [6].) These 
results suggested that the phase of the fcrmion integral must play a crucial role if the SSB 
really takes place. A saddle-point analysis of its effect [7] led to the conclusion that the 
SO(IO) symmetry can be spontaneously broken down to SO{d), where 3 < d < 8. Monte 
Carlo studies have been attempted in Ref. [8] and an intuitive argument for d = 4 has been 
presented. Ref. [9], on the other hand, has provided concrete examples of exactly solvable 
matrix models, in which SO{D) symmetry is spontaneously broken precisely due to the 
phase. 

In this paper, we address the issue of SSB in the IIB matrix model by using the 
Gaussian expansion technique developed in Refs. [10, 11]. Such a method for general su- 
persymmetric models was originally proposed by Kabat-Lifschytz [12] with the particular 
intention of studying the strong-coupling regime of the Matrix Theory [13]. Even at the 

leading order of the expansion, the result turned out to be consistent with the conjectured 
duality to supergravity, and various nonperturbative blackhole dynamics have been ex- 
tracted successfully [14]. However, the possibility of a systematic improvement beyond the 
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leading order was left unclear. Here we will propose a prescription for including higher or- 
der corrections systematically, which indeed yields a rapid convergence to the exact result 
in a simple example. 

We apply this prescription to the IIB matrix model and carry out calculations up to 
the third order. 'Self-consistency equations' allow various solutions which preserve only 
some subgroup of SO(IO). Among them, the solution preserving S0(4) symmetry gives the 
smallest free energy. The value of the free energy comes much closer to the analytic formula 
obtained by Krauth-Nicolai-Staudacher [15] as we increase the order. As a fundamental 
observable which probes the space-time structure, we calculate the extent of space-time in 
each direction. For the S0(4) preserving solution, the extent in the four directions is found 
to be larger than the remaining six directions and moreover the ratio becomes larger as we 
go to higher order. These results provide the first analytical evidence that four-dimensional 
space-time is generated dynamically in the IIB matrix model. 



2. Bosonic Yang-Mills Integral 

In order to illustrate our method, let us consider the bosonic Yang-Mills integral defined 

by 

Z = J dAe-^ , (2.1) 

S = ~Nl3ti[A^,A,]\ ^=^- (^-^^ 

The bosonic matrices A^ {/j, = 1, ■ ■ ■ , D) are NxN hermitian matrices, which we expand as 
Afj, = ^4^ T" with respect to the SU(iV) generators (a = 1, • • • , (iV^ - 1)) normalized as 
tr (r«T^) = ^S"^. The integration measure for A^ is defined as d^ = Y\a=i^ U^=i 
an important consequence of the zero- volume limit, one can actually absorb the parameter 
g by rescaling the dynamical variables Af^^ ^Af^. Therefore, the parameter g is merely 
a scale parameter rather than a coupling constant. The partition function is conjectured 
to be finite [16] ioi N > D/{D — 2), and this conjecture was proved in [17]. A systematic 
1/D expansion has been formulated in [3]. In particular the absence of SO(D) breaking is 
shown to all orders of the 1/D expansion and this conclusion is confirmed by Monte Carlo 
simulations [3] for various D = 3, 4, 6, • • • , 20. The model has been studied by the Gaussian 
expansion assuming that the SO{D) symmetry is not spontaneously broken [10], and the 
numerical results of the VEVs of Polyakov line and Wilson loop [18] have been reproduced 
qualitatively. 

In order to examine the SSB of SO{D) symmetry by means of the Gaussian expansion, 
we repeat the calculation of Ref. [10] without assuming the absence of SO{D) breaking. 
It turns out to be convenient to introduce the rescaled dynamical variables given by 

= P^/^A^ , so that the action takes the canonical form 

S = ~^NtT[X^,X,f . (2.3) 
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The most general SU(Af) invariant Gaussian action without assuming SO{D) symmetry 
can be brought into the form 



^ N 



EI V 
— tr , v^,>0, (2.4) 



by making an appropriate SO(D) transformation. Then we rewrite the partition function 
(2.1) as 

Z = Zo(e-(^-^o))o, (2.5) 
Zo = j dAe-^° = ^-DiN'-m j dXe-^" , (2.6) 

where ( • )o is a VEV with respect to the partition function Zq. Prom (2.5) it follows that 
the free energy F = —InZ can be expanded as 

oo 

F = ^Fk ; Fo = -lnZo, 

k=0 

Fk = -^^{{S-So)'')c,o (forfc>l), (2.7) 

where the suffix 'C in ( • )c,o means that the connected part is taken. The first two terms 
of the expansion are given as 

1 ^ 
Fo = -{N^ - l){l)ln(iV/?V2) - ^In^;^} , (2.8) 

Fi = {S)o - {So)o , (2.9) 
{S)o = l{N'-l)^v,v., (2.10) 

{So)o = liN' - 1)D . (2.11) 

The variational parameters in the Gaussian action (2.4) can be determined in such 
a way that the free energy F calculated up to the first order becomes minimum. This gives 
the self-consistency equations 



Considering > 0, one immediately finds that 



vi = V2 = ■ ■ ■ = VD = y , (2.13) 

which agrees with Rcf. [10]. Thus the Gaussian approximation is able to reproduce correctly 
the absence of the SO(D) symmetry breaking. 
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3. The IIB matrix model 



Let us move on to the IIB matrix model, which is defined by the partition function 

Z = J dAd^e"^ ; S = S^^^ + 5^^) , (3.1) 

where S^^^ is the bosonic action given by (2.2) with D = 10, and 5^^) is given by 

5(F) = -'-Nptr (^'„(r^)«;3[A^, ^p]) . (3.2) 

The fermionic matrices (a = 1, • • • , 16) are traceless N x N hermitian matrices, which 
we expand as ^'q, = T". The integration measure is given by d^' = H^T^ (nl=i d^'"^ . 

The 16 X 16 symmetric matrices are given as = CT^, where are the lOd gamma 
matrices after Weyl projection, and C is the charge conjugation matrix, which satisfies 

= CT^C^ and = C. The partition function (3.1) is conjectured [15] to be finite for 
arbitrary and this conjecture was proved in [19]. 

Again we introduce the rescaled variables X^, and similarly we introduce $q, by $q, = 
0^^^ , so that the fermionic action takes the canonical form 

5(F) = _^jvtr ($«(r^)„;3[^,.,$^]) • (3.3) 

We write down the SU(iV) symmetric Gaussian action without assuming SO(IO) invariance 
as 

So = sfKsf\ (3.4) 
where is the anisotropic bosonic Gaussian action (2.4) and Sq is written as 



where ^ is a 16 x 16 anti-symmetric matrix. 

As before, we introduce Zq, ( • )o and the free energy is expanded as (2.7). Note 
that correlation functions ( • )q q including odd powers of S**^^^ vanish trivially. When we 
evaluate the free energy F by using the formula (2.7), we have to reorganize the expansion 
in such a way that wc can take into account the three-point interaction S**^^^ properly. As 
in Rcf. [11], we consider S'^^\ S^^\ , (S^^'^f to be of the same 'order' and arrive at 
the new expansion 

oo 

F = J2Fk; Fo = -lnZo; (3.6) 



k=0 

k 



^. = -E«w((^(^)-^o)'=-'(s(^))^') , 

' ' c,o 

1=0 



ak,l = k+lCk-i (for k>l). (3.7) 
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This prescription was first proposed in Rcf. [12] for general supersymmetric models ^. 
Tlie zerotli order free energy can be obtained as 



Fo = (iV2 - 1) { 



} , (3.8) 



C-^5^1n^;^-ln(Pf^) 

C = ^{l01n(iV/?^/2) _ i6in(Ar/33/4)| . (3.9) 
The first order correction to the free energy reads 

A = (5(%-(5o)o-^((5(^))%^^. (3.10) 
The first term is the same as (2.10) in the bosonic case. The other terms are given as 

{So)o = -3(iV' - 1) , (3.11) 

{{S^^^f)c,o = l{N^-l)Q, (3.12) 

Q = ^p^v^, p^ = ^Tt[iA-'t^f} , (3.13) 

where the trace Tr is taken with respect to the 16-dimensional spinor indices. Thus, at 
the first order, the free energy is calculated as 



Let us parametrize the 16 x 16 anti-symmetric matrix Aa/S in a SO(IO) covariant way 



as 



^ = ^ 5^ w^,,x B^,x ; ^^.A = c rt , (3.15) 

where w^j^i,\ is a complex totally anti-symmetric rank-three tensor. As it becomes clear 
from this parametrization, the fermionic Gaussian action breaks the full SO(IO) symmetry 
for any nonzero w^^x- The self-consistency equations at the first order read 



= ~Th:{A-^B^,x) + ^^i;^Tv{(^-ir^)2^-ii3^,A} ■ (3.17) 



^ In Ref. [12] the Gaussian action is constructed using a superfield formalism. Similarly we may use the 
four-dimensional superfield formalism regarding the IIB matrix model as the zero-volume limit of 4d A/" = 4 
super Yang-Mills theory. The result turns out to be quite similar to what we obtain below for the solution 
assuming the SO (4) symmetry. 
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In order to go beyond the leading order approximation, we specify the prescription for 
systematic higher order calculations as follows. We choose the parameters in the Gaussian 
action in such a way that the truncated free energy X^^^q is extremized. (Technically, 
we search for solutions by Newton's method near the solutions found at lower orders.) This 
will give the n-th order approximation to the free energy. Similarly to the free energy (3.6), 
the expectation value of an operator O is calculated using the reorganized series expansion 

oo 

(O) = {O)o + Y,0u ; (3.18) 

O.^Y. «M ((^^^^ - Sof-\S^^-^ro)^ ^ . (3.19) 

We truncate the infinite series (3.18) at A; = ra, and evaluate it at the solutions to the ra-th 
order self-consistency equations. 

As a test, we applied this prescription to a toy model [10] with the action S = -^(l>'^- 
Here we can easily proceed up to the 10th order (or even higher) by using Mathematica. 
We found solutions to the self-consistency equation except for the 2nd and 4th orders. At 
each order (other than 2 and 4), we calculated the free energy and the expectation value 
{(j)'^). For the free energy, the discrepancy from the exact result is only 0.1% already at 
order 3, and it becomes 0.004% at order 10. For the expectation value (^^), the discrepancy 
is found to be 0.4% already at order 3, and it becomes 0.03% at order 10. This suggests 
that the above prescription for higher order calculations indeed yields a rapid convergence.^ 
Details shall be reported elsewhere. 

4. Ansatz 

Since the Gaussian action contains too many parameters (10 real numbers from in (2.4) 
and 120 complex numbers from w^i^\ in (3.15)), it seems a formidable task to explore 
the whole solution space of the self-consistency equations. Here we search for solutions 
assuming that SO(d) plus some discrete subgroup of SO(IO) is preserved. For each case 
below (d = 2,4,6,7) the parameters are reduced to two real and one complex numbers. 
For other values of d, we need to keep more independent parameters, and we leave these 
cases for future investigations. 

First let us assume that SO (7) symmetry is preserved, which allows us to set vi = ■ ■ ■ = 
vj = V and all the w^i^x except if8,9,io = w to zero. By further imposing the symmetry 
under cyclic permutations of xg, xg and xio, we set = = viq = v. 

Secondly let us assume SO (6) symmetry, which allows us to set vi = • • • = = 1^ 
and all the w^i^x with indices 1 6 to zero. We further impose the symmetry under 
even permutations of xj, xg, xq and xiq. This requires = vg = vq = viq = v and 

W7,8,9 = -^^7,8,10 = W'T.g.lO = -W'S.Q.IO = W- 

^We have also tried other prescriptions. For instance, when we evaluate the expansions truncated at 
the n-th order, we may use the solution to the first order self-consistency equations. However, the results 
obtained for the free energy and {(/>^) both start oscillating at the third order, and the oscillation becomes 
increasingly violent as we go to higher order. 
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Next let us assume SO (4) symmetry, which aUows us to set vi = ■ ■ ■ = V/j^ = V and 
all the W/ii/A with indices 1 ~ 4 to zero. We also impose the symmetry under a particular 
SO(IO) transformation 

I ^ Xi/ 5 Xi/ I ^ x^ 5 x± I ^ xi , (^"1) 

where (jU, z^) = (5, 6), (7, 8), (9, 10). Furthermore we impose the symmetry under cyclic 
permutations of the three pairs (xs, xq), (xy, xs), {xg,xio). This leads to ^5 = • • • = viq = v 
and w^i^x = w, where /i G {5, 6}, z> € {7, 8}, A G {9, 10}. 

Finally let us assume S0(2) symmetry, which allows us to set vi = V2 = V and all the 
w^uX with indices 1 or 2 to zero. We also impose the symmetry under the transformation 
(4.1) with (/Lt, v) = (3, 4), (5, 6), (7, 8), (9, 10). Furthermore we impose the symmetry under 
even permutations of the four pairs (X3, X4), {x5,xq), (xy, xs), (xg, Xiq). This requires V3 = 
■ ■ ■ = vio = V and w^^x = -Wnup = WfiXp = -WvXp = w, where n G {3,4}, v G {5,6}, 
Ag{7,8}, pG{9,10}. 

5. Results 

For each ansatz preserving the SO(c?) symmetry {d = 2, 4, 6, 7), wc obtain one or two solu- 
tions. In the latter case, we only show results for the one which gives the smaller free energy. 
An analytic formula for the free energy at arbitrary N was obtained in Ref. [15] based on 
explicit numerical evaluations at small combined with other analytical calculations [20]. 
At large AT, the formula gives ^ 

^=l„(v*/) + (ln8-^)+0(!^) . (5.1) 

The Gaussian expansion reproduces the first term correctly for any solution. Therefore, 
we compute the 'free energy density' defined by 

The results are shown in the Table. At the first order, the free energy becomes larger 



d 


/ (order 1) 


/ (order 3) 


p (order 1) 


p (order 3) 


2 


6.49428 


6.50906 


2.17736 


1.49056 


4 


6.15335 


0.74111 


1.85728 


3.37766 


6 


5.75743 


1.54414 


1.87034 


2.24911 


7 


5.52272 


1.62094 


1.95533 


2.15681 



for smaller d. At the second order, we find no solutions to the self-consistency equations. 
This is not so surprising, however, since wc encounter a similar situation with the (f)^ toy 
model as we mentioned below (3.19). At the third order, we find that the free energy 

^The famous factor [20] X]n|jv partition function gives an 0{N~^) contribution in eq. (5.1), 

hence it is irrelevant in the present analysis. 
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becomes minimum at the solution preserving SO (4) symmetry. Note also that the value 
of / obtained for d = 4 comes much closer to the 'exact' result (ln8 — | = 1.32944) as 
we proceed from order 1 to order 3. We also calculate the extent in the /x-th direction 
R/j. = ^ {j^tr Note that Ri = ■ ■ ■ = Rd = R and Rd+i = • • • = Rio = r due to the 

imposed symmetry. At the first order, the ratio p = R/r is given by \fvjv and we find 
that p > 1. At the third order, we observe that the ratio p increases in all the cases except 
for d = 2. 

6. Discussion 

In this paper we have formiilatcd an analytical approach to the spontaneous breakdown 
of SO(IO) symmetry in the IIB matrix model. Our approach is based on the Gaussian 
expansion technique, which was quite successful in the Matrix Theory even at the leading 
order. We have given a prescription for systematic higher order calculations, which, to our 
surprise, yields a rapid convergence in a simple example. 

Here we have carried out our program up to the third order, which can be done with 
reasonable efforts (e.(/., we evaluated 34 four-loop diagrams). We found various solutions 
which break SO(IO) symmetry, and among them the S0(4) preserving solution turned out 
to give the smallest free energy. Moreover the ratio of the extents in the 4 and the other 6 
directions increases as we go to higher order. These results support the conjectured scenario 
that 4d space-time is generated dynamically in nonperturbative superstring theory. Let us 
also recall that the emergence of four-dimensional space-time in the IIB matrix model has 
already been suggested in Ref. [8] based on Monte Carlo results and the branched polymer 
description [2] of its low-energy effective theory. It is encouraging that we arrive at '4d' 
from a totally different approach. 

The analytic formula for the free energy obtained by Krauth-Nicolai-Staudacher [15] 
provides a useful guide-line for the convergence of the present approach. Indeed we observed 
that the free energy calculated for the SO (4) preserving solution comes much closer to the 
KNS result as the order is increased. In the bosonic model, the Gaussian expansion (even 
at the first order) becomes exact in the large-D limit [10]. This fact can be understood 
naturally from the viewpoint of a systematic 1/Z) expansion [3]. In the super symmetric 
models, D is restricted to 4, 6 and 10 and therefore 'the large-D limit' does not make sense. 
Still it is conceivable that the convergence is reasonably fast if D is as large as 10, i.e. the 
case corresponding to the IIB matrix model. The same reasoning explains the success of 
the (leading-order) Gaussian approximation in the Matrix Theory [14]. Nevertheless it is 
certainly worth while to perform higher order calculations and confirm the convergence 
directly. Incidentally we had to restrict the parameter space of the Gaussian action for a 
purely technical reason. It would be interesting to enlarge the parameter space to study 
other SO(c?) preserving solutions, in particular d = 3 and d = 5. 

Finally we would like to comment on the scale parameter g in the action (2.2) and 
(3.2). In the present approach, we rescale the matrices in such a way that the action takes 
the canonical form (2.3), (3.3). This corresponds to setting g^N = 1 in the original action. 
Then the A'^ dependence is eliminated from the self-consistency equations. If (O) converges 
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in high order calculations, it means that the quantity becomes finite in the large-A^ limit 
with g^N fixed. This explains the observed large- A'^ behavior of the space-time extent and 
the Wilson loops in the bosonic and super symmetric models [3, 18]. In the IIB matrix 
model, it is therefore plausible that g^N should be fixed in order to obtain finite Wilson 
loop correlators in the large-iV limit. On the other hand, Monte Carlo results [8] suggest 
that the extent of the space-time (in the 4 directions) may diverge in this large- iV limit. It 
would be interesting to see if such trends appear in the higher order calculations. We hope 
that our analytical approach to the IIB matrix model is useful to extract more information 
on its highly nontrivial dynamics. 
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